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Abstract 

Recently increased accessibility of large-scale digital records enables one to monitor human ac- 
tivities such as the interevent time distributions between two consecutive visits to a web portal by 
a single user, two consecutive emails sent out by a user, two consecutive library loans made by a 
single individual, etc. Interestingly, those distributions exhibit a universal behavior, D(t) ~ , 
where r is the interevent time, and 5 ~ 1 or 3/2. The universal behaviors have been modeled via 
the waiting-time distribution of a task in the queue operating based on priority; the waiting time 
follows a power law distribution P^{t) ^ with either a = 1 or 3/2 depending on the detail of 
queuing dynamics. In these models, the number of incoming tasks in a unit time interval has been 
assumed to follow a Poisson-type distribution. For an email system, however, the number of emails 
delivered to a mail box in a unit time we measured follows a powerlaw distribution with general 
exponent 7. For this case, we obtain analytically the exponent a, which is not necessarily 1 or 
3/2 and takes nonuniversal values depending on 7. We develop the generating function formalism 
to obtain the exponent a, which is distinct from the continuous time approximation used in the 
previous studies. 



PACS numbers: 89.75.Hc, 89.70.-a, 89.20.Ff 



I. INTRODUCTION 



In the digital era, human activities can be easily monitored and quantified by analyzing 
digital records such as the dates of sending or replying to emails, and financial transactions. 
Interestingly, human activities generate emerging patterns: the interevent time distribution 
of human activities follows a power law, and its exponent is either 1 or 3/2 in many cases [l|,y, 
y, y, [sj. Such a bursty nature of human dynamics has been understood to be a consequence 
of queuing processes driven by human decision making. Barabasi introduced a queuing 
model operating in the priority-based protocol At each time step, a task arrives at such 
a queue and is assigned a priority Xi chosen randomly from a distribution p(x). Then, with 
probability p, the task with the highest priority is selected for execution and removed from 
the list. With probability 1 — p, a task is randomly selected irrespective of its priority and is 
executed. This model was successful in analytically reproducing the empirical result 
the waiting time of a task in the queue before being executed, which is denoted by r, follows 
a power-law distribution Pw(t) ~ r^^. The result is independent of distribution p{x). The 
power law Pw(t) ~ r"^/^ is reproduced by allowing the queue length to vary in time 

To analyze both fixed-length and fiexible-length queues, the Barabasi's model was ex- 
tended as follows. In each time step, a task arrives with probability A, and the task with 
the highest priority in the queue list is executed with probability fi. Operation of this queue 
system is schematically shown in Fig. [TJ^a). Since the dynamics of the queue is stochastic 
ifO<A<lorO</i<l, the queue length generally changes in time. This model is a 
type of the M/G/1 queuing system with a priority selection rule proposed in the seminal 
work of Cobham in 1954 . This model was analytically studied recently. The waiting-time 
distribution of a task in the queue changes depending on A and fi. (i) When A = /i = 1, the 
number of tasks in the queue is fixed, and the waiting time of tasks obeys Pw(t) ~ 3]- 
(ii) When A = p < 1, P„(r) ~ r-3/2 (iii) When A < /i < 1, Pw(r) ~ r-3/2g-r/ro 
for T <^ tq and Pw(t) ~ r'^/^e"^/"^" for r ^ tq, where the characteristic time scales as 
To = 1/ (y/fi- — [8|. (iv) When /i < A < 1, tasks with priority x < (A — /i)/A wait in the 
queue forever without being executed. Tasks with priority x > (A — /i)/A are executed with 
the waiting time r following Pw(t) ~ r^^/'^ 8|]. 

Previous studies focused on the case in which incoming tasks are independent of each 
other and delivered to the queue at a constant rate. Thus, the number of incoming tasks in 
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a unit time follows the Poisson distribution. This is the case observed in, for example, the 
number of requests for wireless phone calls arriving at a cell station in a unit time (see inset 
of Fig. [2]). However, we observe that the number of emails received by a single user in a 
unit time is heterogeneous and follows a power-law distribution (see Fig. [2]) . Time intervals 
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between consecutive tasks arriving at a server computer 

user's hypertext markup language (HTML) requests js], which are closely related to the 
number of incoming tasks per unit time, also show similar patterns. The origin of such non- 
uniform numbers of incoming tasks is not known yet, but may be consequences of multiple 
correspondences with multiple people or self-similar patterns in the number of data packets 



arriving at a given router 



141 ]. Such bursty arrivals of tasks may significantly change the 



behavior of priority queue systems. For example, a more skewed distribution of the number 
of incoming tasks per unit time may result in a more skewed waiting-time distribution of a 
task Pw(t), as briefly suggested in [12]. In this paper, we study the waiting-time distribution 
of a task in the queue for the case of heterogeneous numbers of incoming tasks. We find 
that the universal power- law exponent a = 3/2 for P„(r) ~ r^" occurs as a limited case 
and obtain other values of a depending on the power-law exponent of the distribution of the 
number of incoming tasks. 



II. MODEL 



We study the queue model defined as follows: in each discrete time step, n tasks are 
delivered to the queue, where n is distributed according to a power law A„ = Xn~'^/({'j) 
{n > 0), Ao = 1 — A, where < A < 1 and ^(7) = X]^=i '^'"''^ Riemann C, function. 

Each task is assigned a priority x uniformly distributed on [0,1]. At the same time, the task 
with the highest priority in a queue is executed with probability /i (0 < /i < 1). Operation 
of this queue system is schematically depicted in Figs. [I](b) and[T](c). The queue length 
is unbounded so that the queue accommodates all incoming tasks. This model generalizes 
the model introduced by Grinstein and Linsker (GL) [sl, which corresponds to Aq = 1 — A, 
Ai = A, and A.„ = for n > 2 in our model. 

We will obtain the waiting-time distribution P„(r) for a task in the queue. To this end, 
we start with the probability that there are m tasks with priority larger than or equal to x in 
the queue at time t, which is denoted by Qx{rn,t). We denote the queue-length distribution 
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FIG. 1: (Color online) Schematic representation of queueing protocols, (a) A queue system pro- 
posed by Grinstein-Linsker, in which at most one task (filled circle) arrives in the system per time 
step, (b) The queue system we consider in this paper, in which input tasks (filled circles) can be 
bursty, (c) Operation of the queue system shown in (b): At time to, there are m tasks (black 
circles) with priority > 2; in the queue. At time t^ + l, n tasks (black and gray circles) arrive in the 
queue with probability A„. Among them, n — i tasks (black, not gray, circles) have priority > x. 
This event occurs with probability (^'^J (1 — x)^~^ x*, where < i < n. The task with the largest 
priority is executed with probability ^. No task is executed with probability 1 — /i. At time to + t, 
the queue does not contain any tasks with priority > x for the first time. 

in the steady state by Qxim) = limt^oo Qx{i^,t). Note that the steady state exists only 
under a certain condition, as discussed later. We define Gx(m, r) to be the probability 
that a given task with priority x arriving in the queue at time t = to is executed at time 
t = tQ + r. When the task arrives in the steady state, there are already m tasks in the 
queue with priority larger than or equal to x, where m is distributed according to Qx{fn). 
All of these m tasks are executed before the given task is executed. Then, the waiting-time 
distribution is obtained via the following formula 81]: 

00 „i ^ 

X 

{m,T), (1) 

where Gx{m,T) is equivalent to the first passage probability that a random walker starting 
from position m > arrives at the origin at time r for the first time. For a constant rate 
of incoming tasks, Qx{m), Gx{m,T), and P„(r) can be obtained explicitly j^. However, 
due to the complexity of our problem, we obtain them implicitly in terms of the generating 
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FIG. 2: (Color online) Distributions of the number of incoming tasks. The main panel shows the 
number of tasks delivered to an email box of an anonymous user per unit time [9|, which follows 
a power-law distribution with slope —1.5. Different lines correspond to different bin sizes, namely, 
500 (O)) 800 (A), and 1000 (□) seconds. Note that the slope —1.5 is not universal. It depends on 
users and can be as small as —3. We chose a user with the largest dataset. (Inset) The number 
of wireless phone calls arriving at a cell station in 10 seconds for the peak time (i.e., 12:00-20:00) 
(A) and for the entire day (O)- Both data fit well to the Poisson distribution (black solid line), 
which decays even faster than the power law with exponent —6 (dotted line). 



functions. We define the generating function 

oo 

P^(s)^^P„(r)s^ 

where < s < 1. Then, 



T=l 



OO „1 ^ 

^w(s) = ^ / dxQ^{m)g^{m,s) 



(2) 



(3) 



where Qxijn.s) = ^^^(m, r)s'^. Because the number of tasks in the queue decreases at 
most one per unit time, we obtain 

G,{m, t) = Y, G.{m - l,t - t)Ut), (4) 

T 

where fx(t) = Gx{l,t). Equation (j4]) is expressed in terms of generating functions as 



Qxim, s) = Qx{rn - 1, s)J^^(s) 



(5) 
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where J-'x{s) = fx(t)s^. Applying Eq. ([5]) repeatedly, we obtain 



Then, Eq. ([3]) is written as 

oo „i ^ 

m=0 



(6) 



(7) 



where Qx{z) = Y.m.=oQxi^)^"'- 

Once we derive Qx{z) and J-'x{s) explicitly, we obtain the waiting-time distribution of a 
task in the queue, namely, P^^t). We will show that the waiting time exhibits a power-law 
behavior Pw(t) ~ t~", where the values of a are shown in Table [B The analytic solutions 
are confirmed numerically in Fig. [3l Using our generating function formalism, we can also 



reproduce the results derived in Ref. 



as shown in Appendix IVI[ 
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FIG. 3: (Color online) The waiting-time distribution Pwit). (a) The case {n)x < f-t. Given A = 0.3 
and /_f = 1.0, shown are numerically obtained PwiT) for 7 = 2.5 (□), 3.0 (O)) 3.5 (A), and 4.0 (0), 
yielding to {n)x 0.58, 0.41, 0.36, and 0.33, respectively. Solid lines indicate Pw{t) ~ t~^'^~^\ (b) 
The case {n)x > /-f with 2 < 7 < 3. Given A = 0.5 and fi = 0.5, shown are numerically obtained 
P„(t) for 7 =2.1 (□), 2.5 (Q), 2.8 (A), and 3.0 (0), yielding {n)x ~ 3.39, 0.97, 0.75, and 0.68, 
respectively. Solid lines indicate PwiT) ~ r~^'^"'^'^^^^'^~^\ (c) The case {n)x > n with 7 > 3. Given 
A = 0.5 and ^ = 0.3, shown are numerically obtained Pw{t) for 7 = 3.3 (□), 3.8 (O)) 4.0 (A), and 
4.5 (0), yielding {n)x ~ 0.62,0.57, 0.56, and 0.53, respectively. The dotted line is a guideline with 
slope —1.4, close to the theoretical value —1.5. 
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TABLE I: Power-law exponent a of the waiting-time distribution P^{t) ~ t 

Wa < ("-)a > M 

2<7<3 7-1 ^ 
7>3 1-^ i 

III. THE QUEUE-LENGTH DISTRIBUTION 

In this section, we calculate the queue-length distribution in the steady state by using 
the generating function Qxi,^^- The master equation for Qx{rn,t^ is given by 

oo 



g^(m,t + l)=^^AjX^Q^(m + l,t) (8) 

m oo ^ .s 

+5^ (1 - /i) 5^ A, ) (1 - xf x^-'QM - 2, t) (9) 

m oo / ■ \ 

+E/^ E A,(^;^ J(l-xr+^x^--U(rn-z,t) (10) 

m 

= Pm-i^mQx {m-i,t), (m > 1). (11) 

i=-l 

In the above equation, the three terms in the right-hand side (RHS) correspond to different 
types of events that occur in a unit time. The first term ([8]) represents the case in which j 
(j = 0, 1, ■ ■ ■) tasks arrive in the queue with probability Xj, the priorities of all j tasks are 
smaller than x, and one task is executed with probability fi. The second term (Q represents 
the case in which j (j = 0, 1, ■ ■ ■ ) tasks arrive in the queue with probability Aj, i tasks out of 
the j tasks have priorities larger than or equal to x, and no task is executed with probability 
1 — fi. The third term fllUI) represents the case in which j (j = 0, 1, . . .) tasks arrive in the 
queue with probability Aj, i + 1 tasks out of the j tasks have priorities larger than or equal 
to X, and one task is executed with probability fi. For later discussion, we denote by Pm-i^m 
in Eq. (fTTj) the transition probability of the random walk from position m — i to position m 



7 



in a unit time. The master equation at the boundary is given by 

oo 

g,.(o,t + i)=/i5^A,x^g,(i,t) 



j=Q 



+ 



j=0 j=l 



X) X-' 



X-' 



j=0 



Qx(o,t) 



(12) 



Based on Eqs. (!8|)- (fT2|) . we calculate the generating function Qx{z) for the steady-state 
queue-length distribution Qxijn) = lim^^oo Qx{^, t). Specifically, the generating function of 
Eq. ([H]) is equal to /iA(x) Qx{z) — Qx{0) jz in the steady state, where k{z) = Xljlo -^i-^"'- 
The generating function of Eq. ([9]) is equal to (1 — ii)Qx{z)k [(1 — x) z + x\. The generating 
function of Eq. ffTOj) is equal to fiQ.x{z) {A[{1 — x) z + x] — A (a;)}. The generating function 
of Qx{0) in Eq- (fT2ll is equal to fiA{x)Qx{0). Combining all these terms, we obtain 

/i4(0)(z-l)A(a;) 



^"'^^'^ z-{n + z- fiz)A [(1 -x)z + x] 

To eliminate Qx{0) from Eq. ( fT3l) . we exploit the condition 
denominator and the numerator of Eq. f[T^ converge to zero as z 
L'Hospital rule to Eq. (IT^ to derive 

Qx{0) = [fi-{l-x) {n),]/{fiA (x)), 

where {n)x = Xl^o""^"- Plugging Eq. f|T^ into Eq. f|T3|) yields 

[/i-(n)A(l-a;)](z-l) 



(13) 

: 1. However, both the 
1. Thus, we apply the 

(14) 
(15) 



z — {fi + z — fiz)A [(1 — x) z + x] 
For the steady state to exist, the incoming rate of the task with larger than or equal to x 
(i.e., {n)x{l —x)) must be smaller than the execution rate n [l, 8]; Ai = fi— {n)x{l —x)>0 
is required. In addition, {n)x must be finite, which is equivalent to the condition 7 > 2. 
The mean queue length denoted by {m{x))Q is derived as 

dQx(z)' 



dz 



z=l 



2(1 - fi){n)x{l -x) + {{n^)x - {n)x){l - xf 
2Ai 



(16) 



Equation (fT6|) implies that {m{x))Q diverges when {n'^)x does, that is, when 7 < 3. When 
7 > 3, the queue length is finite for x = if and only if > {n)x and diverges as 1/ (/x — {n)x) 
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as {n)x approaches /i from below, which extends the results in 8|. As x — and {n)x fi, 
(m{x))Q diverges as 1/x, which is also consistent with the previous result Q]. 

To calculate the asymptotic behavior of the steady-state queue- length distribution Qxijn), 
we assume /i > {n)x{^ — x) and 7 > 2, for which the steady state exists. When 2 < 7 < 3, 



A.{z) is expanded near 2; ^ 1 as follows 15|: 



A(^) = 1 - {n)xil -z)+ c,(l - zy~' + o ((1 - zy~') , (17) 
where is a constant. Inserting Eq. f|T7j) into Eq. f|T5|l leads to 



For 3 < 7 < 4, we obtain 

A(^) = 1 - {n)x{l -z) + _ ,f _ ^^^^ _ ^y-i + o ((1 - zy^') , (19) 
which leads to 

Qx{z) = 1 + (m (x))q (. - 1) + _ ,)7-2 + ^ ((1 _ ^)7-2) (20) 



Similar expansions hold true for 7 > 4. By applying the Tauberian theorem [15| to Eqs. (fT8|) 
and (120|) . we obtain 

4(m)~^ (m^oo) (21) 

for 7 > 2. Equation fl^Tl) is consistent with the result under the first-in-first-out (FIFO) 
protocol |l6|. This is because, when fi ^ ('^)a(1 — x), tasks are executed upon its arrival in 
the steady state so that the priority-based protocol can be regarded as the FIFO-based one. 



IV. FIRST-PASSAGE PROBABILITY 

In this section, we derive J-'x{s) = Ylt^i fx(t)s^ = Ylt^i^x{^,t)s^- Recall that Gxijn.t) 
is the probability that a given task with priority x is executed at time t after its arrival, 
provided that there are m tasks in the queue with priority larger than or equal to x when 
this task arrives. This quantity can be interpreted as the first passage probability that a 
random walker on a half line starts from position m and arrives at the origin at time t for 
the first time. The probability that the random walker moves from i to j in a unit time is 
given by Pi^j [see Eq. (fTTj) ]. 
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The generator of the one-step transition of the random walk before reaching the origin is 
represented by 

oo 

V{Z) = ^ Pm^m+iZ' = (l-fI+^^A[{l-x)z + x]. (22) 
i=-l 

Note that the RHS of Eq. (122|) is independent of m because the transition probabihty is 
homogeneous in space. 

The amount of a single jump that the random walker makes to the right is unbounded, 
because it is equal to the number of incoming tasks with priority larger than or equal to x. 
However, the amount of a jump to the left is at most one, which yields a useful relation, 

G.{hs)=J',{sy. (23) 

Using Eqs. fl22l) . fl23|) . and the recursion relation 
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18| 



Ut) = pi^o + Pi^iG^il, t-l)+ Pi^2G,{2, t - 1) + . . . , (24) 
we obtain the following self-consistent equation for the generating function: 

oo 



1=0 
oo 



= s^pi^ij^^{sy 

= s[(l-/i)^,.(s)+/i]A[(l-x)^,(s)+x]. (25) 

The first s in the RHS comes from the unit time spent by a single transition starting from 
m = 1. After this transition, the generating function of the number of tasks with priority 
larger than or equal to x in the queue is zV{z). Since each such task incurs an execution 
time distributed according to {fxif)}, we replace z of zV{z) by to obtain Eq. (!25l) . 

We evaluate J^x{s) in the limit s ^ 1 using Eq. fl25l) . To guarantee that the task with 
priority x is eventually executed, = 1) = 1 has to be satisfied. To check if this condition 
is fulfilled, we put s = 1 in Eq. fl25l) to obtain 

T^il) = s [(1 - /i) (1) + /i] A [(1 - x) Tx (1) + x\ . (26) 

The left-hand side (LHS) and the RHS of Eq. (l26ll are plotted in Fig. |4]as functions of J^z(l), 
where J-'x{^) is regarded as a variable for the sake of this analysis. Note that the RHS of 
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Eq. (126|) is positive at = 0. Figure H] implies that Eq. (126|) has the unique solution 

JF^(I) = 1 if and only if the slope of the RHS of Eq. fl26l) at J-'x{^) = 1 is less than or equal 
to unity, that is, 

d 



~ ^^^^^ ^^^^ " + 1^41)=! - 

Equation ( 1271) is equivalent to Ai > 0, which is what we already assumed. 



(27) 




1 ^.(1) 



FIG. 4: (Color online) Schematic representation of the LHS and the RHS of Eq. (I26p as functions 
of .7^,(1). 



In the following, we obtain the solution of the self-consistent equation ( l25l) by assuming 
that fx{t) follows a power law. 

Case (i): fi > {n)x- In this case, Ai > holds for all x. When 2 < 7 < 3, combining 
Eqs. (HZD and ([25]) yields 



1 + - 1) + ^^^^^^4^(1 - sy-' + o ((1 - sr') . 



(28) 



When 3 < 7 < 4, combining Eqs. (fT9l) and (125!) yields 



c^(l - x)T-i 



l-s)^-^ + o((l 



\7-l 



)(?9) 



where = ((7^^)a - ('^)a)(1 - xy/2 + (1 - iJ,){n)x{l - x) > 0. Note that the coefficient of 
(1 — s)2 is positive. In a similar manner, we can show for 7 > 4 that the leading singular term 
of J-'xis) is equal to (— l)'^'^^~^c^(l — x)^~^{l — s)^^^/Al, where [7] = min{i; i > 7, i G Z}. 



Thus, we obtain fx{t) r-^ t ^ with /5 = 7 for 7 > 2 using the Tauberian theorem [15 1 
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Case (ii): /i = {n)x. Because Ai = for x = 0, we cannot apply the results obtained for case 
(i). For example, l/Ai = 1/ {{n)xx) in the coefficient of (s — 1) in Eq. fl28|) diverges as a; ^ 0, 
implying that the exponent f] is smaller than 2 near x = 0. Actually the long-time behavior 
of fx{t) is dominated by the tasks whose priority is near x = Q]. Thus, we assume 

Txis) = 1 - c^(l - sf-' + o((l - sf-') (30) 

with 1 < (3<2. 

When 2 < 7 < 3, the RHS of Eq. ([23]) is written as 

= s [(1 - /i) j;(s) + /i] A [(1 - x)J^x{s) + x] 
= sTx{s) + /iS [1 - J^x{s)] + 

s [ {n)x (1 - x) [J-.(s) - 1] + c,(l - xf-' [1 - ^.(s)](^-^) + ■ ■ ■ } (31) 

Plugging Eq. ([30D into the LHS and RHS of Eq. ([25]) leads to 

(1-.) [l-c,(l-s)^-^ + ---] 
= {n)xxcp{l - sf-^ + c^(l - xy-^cl-\l - s)('3-^)(^-^) + ■ ■ • . (32) 

If {n)xx > (1 - s)('>'-2)/{7-i)^ the first term of the RHS of Eq. ([32]) is much larger than 
the second term as s 1 so that [3 = 2 and cp = 1/ {{n)xx). Conversely, if {n)xx <^ 
(1 — sY'^~'^^^^"'~^\ the second term dominates the first term so that /3 = 1 + 1/(7 — 1) and 

C^ = C-'^^'-'^/{l-x)c:.C-'/'~"'-'\ 

When 3 < 7 < 4, as in the case of 2 < 7 < 3, Eqs. ([25]), and ([30]), with an 

appropriate assumption of 1 < /5 < 2, yield 

{l-s)+ 0(1 -s) = {n)xxcp (1 - sf-' + A2C%1 - sf^^-'^ 

- c^(l - xy~\}-\l - s)('3-i)(7-i) + . . . . (33) 

If {n)xx ^ ^2(1 — s), the first term in the RHS of Eq. (133]) is much larger than the second 
term. Then (3 = 2 and cp = 1/ {in)xx). Conversely, if {n)xx <ti 1/^2(1 — s), the second 
term is much larger than the first term so that /3 = 3/2 and Cjs = l/^fA^. The third term is 
always much smaller than the second term as s ^ 1. 

Case (iii): /i < {n)x- The task in the queue accumulates at rate {n)x — /i. In this case, only 
the tasks with priority x > xm = {{n)x — jj) /{n)x are executed, and the analysis can be 
ascribed to case (ii) |8|. Distributions of the priority of tasks in the queue in the steady 
state are shown in Fig. [5] for some values of {n)x and [i. 
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FIG. 5: (Color online) Distributions of the priority of tasks in the queue in the steady state. 
Distributions of x for different sets of {n)x and /i are shown. Only the tasks with priority x > 
xm = {{n)x — fi)/{n)\ are executed, (a) A = 0.5 and /i = 0.5. We set 7 = 2.1, 2.5, 2.8, and 3.0, 
which yield xm = 0.85, 0.48, 0.33, and 0.26, respectively, (b) A = 0.5 and ^ = 0.3. We set 7 = 3.3, 
3.8, 4.0, and 4.5, which yield 0.52, 0.47, 0.46, and 0.43, respectively. In each panel, the four plots 
almost collapse onto one. 

V. THE WAITING-TIME DISTRIBUTION 



Using Eqs. ([7j), (ITSl) . and J-'x{s) we obtained for the three cases, we calculate the waiting- 
time distribution as follows: 

Case (i): /i > {n)\. The leading singular term of Qa;(-^a;(s)) is equal to (—l)'^'^^c^(l— a;)'^~^(l — 
s)'^~^/yl5^~^. Then, we obtain 

V^{s) ^ {I - sV-\ (34) 
which yields Pw(t) ~ t'^'^'^'^ for 7 > 2. 

Case (ii): fi = {n)x- In this case, we use Eq. fl3Ql) with values of /5 and depending on 7 
and X. 

For 2 < 7 < 3, we obtain 



Vjs) 



7-2 

dx{n)\xc^~' (1 — s)t-i + / + 

i(l-s)^ 



Therefore, Pw{t) ~ t 



WA_Z 1 I (l_s)(.-i) + 



^(27-3)/{7-l) 



(35) 
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For 3 < 7 < 4, we obtain 



VM^\ dx^M^=+l dx 



1 - S) 

Therefore, P„(r) ~ r~^/^. Similar calculations yield P„(r) ~ r^^^'^ for 7 > 4. 

Case (iii): /i < {n)x- Since the analysis can be ascribed to case (ii), we obtain -Pw(t) ~ 

^-(27-3)7(7-1) for 2 < 7 < 3 and P„(r) ~ r"3/2 for 7 > 3. 



VI. DISCUSSION AND SUMMARY 

The analytic results are summarized in Table [T] and confirmed numerically in Fig. [31 The 
power-law behavior of the waiting-time distribution -Pw(t) ~ can be diverse in that a 
can take general values, rather than a = 1 or 3/2. Consistent with this, the intertransaction 
time of a stock broker obeys the power-law distribution with a ~ 1.3 with an exponential 
cutoff [31]. 

Our results are compatible with those derived from the continuous time approximation 8] 
and the fractional derivative 19|. The generating function approach that we have developed 
can be useful for studying further problems. For example, we show in the Appendix that 
our approach considered in the limit A, /i — > reproduces the results for the GL model Q]. 
Furthermore, GL as well as we are successful in deriving the exponential cutoff for A < 
as Pw{t) ~ r'^/^e"'^/'^'' with tq = l/{y/JI — V^)^. However, for the model with general 
distributions of the number of incoming tasks, the explicit form of the exponential correction 
factor is not obvious. 

In our priority queue model, the jump distance of the equivalent random walk is 
unbounded to the right, whereas it is at most one to the left. In real queue systems, 
however, more than one tasks may be executed in a unit time. Therefore, a natural 
extension of our model is to allow the number of executed tasks in a unit time to exceed 
one. To be specific, in addition to the heterogeneity of the number of incoming tasks, i.e., n 
tasks are incoming with probability A„ ~ n~'>''" per unit time, we can suppose that i tasks 
are executed with probability fie ~ Our numerical results for P„{t) seem to fit the 

formulas shown in Table [H with the exponent 7 replaced by the minimum of 7in and 7out , 
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as far as both and 7out are larger than 2 (not shown). This suggests that the dominant 
tail determines the behavior of the waiting-time distribution in the priority queue system. 
In particular, when the distribution of the number of executed tasks is neither binary nor 
heavy-tailed (e.g. purely exponential), which may be true for many real queues, our results 
hold because 7 = 7in- 

Acknowledgments This work was supported by the KOSEF grant for Acceleration Research 
(CNRC) (Grant No. R17-2007-073-01001-0), the KRCF, and Grants-in-Aid for Scientific 
Research from MEXT, Japan (Grants No. 20760258 and No. 20540382). 



COMPARISON OF THE GRINSTEIN-LINSKER SOLUTION AND THE 
GENERATING-FUNCTION SOLUTION 

GL analyzed a priority queue model in which, in a unit time, a new task arrives with 
probability A and the task with the highest priority in the queue is executed with probability 
/i, which corresponds to Aq = 1 — A, Ai = A, and A„ = for n > 2 in our model |8|. 
They obtained the solution of the waiting-time distribution by analyzing the continuous- 
time dynamics. We compare the GL solution and the solution derived via the generating 
function in the GL limit. 



1. The queue-length distribution 

The generating function of the queue-length distribution in the steady state is given in 
Eq. ( IT5|) in the main text. By substituting A{z) = 1 — \ + \z and {n)^ = A into Eq. (ITSll . 
we obtain 

Q (z) - f^-Kl-x) 



which leads to 



~ /i-A(l-x) f \{1 - - x) \ 

= (T^aTa^)^ ItT^aTW ' ■ ^ ^ 



Using the continuous-time approach, GL derived 



= f^il-£) | . (39) 
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Equations (138|) and (139|) are consistent in the limit A,/i — > 0. 



2. The waiting-time distribution and the exponential cutofF 



To obtain the waiting-time distribution of a task, we use the following theorem [20|, l2l|, |22| : 



Theorem: Suppose that, for real numbers s* and JF*, a power series J^{s) = Xlt^i 0'{t)s'^ 
with nonnegative coefficients a(l), a(2), . . . satisfies the following equations fHOl) . fHT]) . and 

is analytic near = (40) 

if |s| < sM^I < J^*, F(s,^) = ^^^^^ = Oif andonly if (s,^) = (s*,^*); (41) 

^^^0, q^^O, (42) 

Then, 

To apply this theorem to the GL queue model, we define 

F (s, ^) = s [(1 - /i) + /i] {1 - A + A [(1 - x) + x] } - J^, (44) 

so that F (s, J^x (s)) = holds, where J-'x{s) = J2t^i fx(t)s^ is the generating function of the 
first-passage time probability. Then the other main condition of the theorem [see Eq. fHTl) ] 
reads 

^^^^ = s(l - /i) {1 - A + A [(1 - x) ^ + x]} + s [(1 - /i) ^ + /i] A(l - x) - 1 = 0. (45) 
The solution to Eqs. (jHj) and (H5!) with the minimum absolute values is given by 

^ (46) 



1 - A - /i + 2A/i + Ax - 2A/ix + 2^A(1 - A + Ax)(l - x)/i(l - fi) ' 
^ / /i(l - A + Xx) . . 



The rest of the conditions of the theorem are satisfied with s* and JF* given by Eqs. (14611 
and (1471) . Equation (l43l) implies that the tail of the first-passage time probability decays as 
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m « t-^'^s-'. This asymptotic is also derivedbvditectly calculating - (s' - sf'^ 



as s t s* and using the Tauberian theorem 15|, l23|, l24l |. 



The generating function of the waiting-time distribution of a task is equal to that of 
the queue-length distribution given by Eq. f|T5|) with s replaced by To calculate 

the asymptotic of the waiting-time distribution, we erase A by combining Eq. (fT5|) with s 
replaced by J-'x{s) and Eq. ( fTSl) . which yields 

Inserting Eq. (148|) into Eq. (l44l) results in 
(1 - A + Ax)/i(l - s)Ql {Tx (s)) + [(1 - A + Ax + /i) s - 1] Q,. (J", (s)) - A\s = 0. (49) 



Applying the theorem (H0l) - (H3l) with T = [Tx (■§)) leads to the same equation (H6l) . 
Therefore, the waiting-time distribution has the same asymptotic as the first-passage time 
probability, that is, Pw{t) ~ r~^/^s*~^. This asymptotic is also derived by solving Eq. fH9|) 
as V^{s) ~ (s* - as s t s*. 

To evaluate s*, we denote the denominator of the RHS of Eq. (HGl) by H{x), with A and 
fixed. The existence of the exponential cutoff in the first-passage time and the waiting-time 
distribution is equivalent to H{x) < 1 (0 < Vx < 1). 

A straightforward calculation yields (PH/dx^ < 0, Yrnvx^i dH / dx = —oo, and that 
dH/dx = has a unique solution x = {X — fi) / X. As explained in the main text and 
in previous literature j8|, the analysis of case A > yU is ascribed to that of case A = /i. 
Therefore, we assume X < fi and obtain dH/dx < (0 < x < 1). Then the maximum of 
H{x) is realized at x = 0, so that the smallest s* is equal to 

= 77^ = ^ , . (50) 

H{0) 1- X- 12 + 2X12 + 2^X{1 - A)^(l - /i) 

When /i = A, we obtain s* = 1. The asymptotic of the waiting-time distribution is 
-Pw(t) ~ r~'^/^, which is consistent with the results in js] and coincides with our results for 
7 > 3. 

When /i > A, we obtain s* > 1 and -Pw(t) ~ -j--3/2^-t/to ^ where Tq = 1/lns*. In the 
limit A, /i ^ 0, our discrete-time model tends to GL's continuous-time queue dynamics. By 
inserting A = A'Ar, fi = /i'Ar, and r = t'/At into Eq. ( l50i) and letting At — * 0, we obtain 
Pw{t) ~ j-'-'ifi^-r'/To^ where Tq = l/(v^ ~ VX)"^- The predicted Tq agrees with the one 
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derived by GL. They concluded P„(r) ~ 7-'--''/2g-r'/ro r < tq and P„(r) ~ r'-^/^g-r'/ro 
for r ^ tq. Our results only reproduce the asymptotic on the intermediate timescale (i.e., 
r ^ To) because tq diverges as A, — 0. 
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